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Abstract 

^ ' Based on the real-time formalism, especially, on Thermo Field Dynamics, we 

■ derive the Schwinger-Dyson gap equation for the fermion propagator in QED and 

Four-Fermion model at finite-temperature and -density. We discuss some advantage 
of the real-time formalism in solving the self-consistent gap equation, in comparison 
with the ordinary imaginary-time formalism. Once we specify the vertex function, 
we can write down the SD equation with only continuous variables without per- 
forming the discrete sum over Matsubara frequencies which cannot be performed 
in advance without further approximation in the imaginary-time formalism. By 
solving the SD equation obtained in this way, we find the chiral-symmetry restor- 
ing transition at finite-temperature and present the associated phase diagram of 
strong coupling QED. In solving the SD equation, we consider two approximations: 
instantaneous-exchange and po-independent ones. The former has a direct corre- 
spondence in the imaginary time formalism, while the latter is a new approximation 
beyond the former, since the latter is able to incorporate new thermal effects which 
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has been overlooked in the ordinary imaginary-time solution. However both approx- 
imations are shown to give qualitatively the same results on the finite-temperature 
phase transition. 
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1 Introduction 



Quantum chromodynamics (QCD) exhibits very characteristic dynamics of confinement 
of quarks and gluons in long distance as well as the asymptotic freedom at short distance. 
According to the standard cosmology, however, our universe ought to have undergone the 
confinement-deconfinement transition of quarks and gluons in the early stage of evolution. 
This implies that the hadron phase at present should transfer into the quark-gluon-plasma 
phase, if the universe goes back to the past at high-temperature and -density m, ||, ||, Q 
Furthermore, the quark-gluon-plasma is not merely a byproduct of the theory in the sense 
that its existence will be verified by experiments of heavy ion collision in the following 
several years 

At high temperature, the effective coupling constant of QCD gets small owing to the 
asymptotic freedom p. Hence it is natural to consider that the perturbation theory can 
be applied to QCD at high temperature successfully @, ^ ||, 0]. However it becomes 
evident in the early stage of the investigation of high-temperature QCD that infrared 
(IR) singularity prevents one from calculating higher order terms in perturbation theory 

In view of this, our study aims at performing the non-perturbative analysis on the 
phase transition at finite-temperature and -density (FTD). For this we study the self- 
consistent (SC) equation represented by the Schwinger-Dyson (SD) equation and the 
Bethe-Salpeter (BS) equation. Such SC equations are simultaneous non-linear integral 
equations and are able to incorporate various non-perturbative effects self-consistently 
1^. The concept of the SC equation is not new and indeed the SC equation is a well- 



known traditional method to solve the many body problem in condensed matter physics 
and nuclear physics [|T6|, [l^. Nevertheless SC equations turn out to provide the powerful 
method also in elementary particle physics: Recently various types of SC or gap equations 
have been extensively studied to investigate the dynamical symmetry breaking in gauge 
field theories and have succeeded to reveal rich phase structure of strong coupling gauge 
theories |TB[ . In this paper we extend this type of research to the case of finite-temperature 



and density. 

We desire to develop the method to study the phase transition at finite-temperature 
and -density from the quantum field theoretical point of view beyond the thermodynamic 
or phenomenological treatment. For this, we need at first the formulation of FTD field 
theory. Traditional Imaginary- Time Formalism (ITF) due to Matsubara [0 is very effi- 



cient to develop the perturbation theory of FTD field theory. In ITF, the field theory at 
zero-temperature can be transferred into the FTD field theory according to the following 
procedures: in a given Feynman diagram, each integral over the temporal component 
of a "fermion" loop momentum is replaced by an infinite sum over "odd" Matsubara 
frequencies according to the prescription: 

/ v^/(Po) fiPo - (2n + l)nT - ^/x), (1) 

where n denotes the integer, T the temperature and fi the chemical potential. In contrast, 
each integral over the temporal component of "bosonic" loop momentum is evaluated by 
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summing over "even" Matsubara frequencies, 



/ -^g{k,) Y: 9{h - ^n-nT). (2) 

However applying this method to the non-perturbative study is sometimes either conve- 
nient nor easy. Indeed the infinite sum is very cumbersome particularly in solving the SC 
equation, because the convergence of the infinite series is problematical and the discrete 
sum cannot be performed in advance before we know the solution, in sharp contrast with 
the perturbation theory where all the propagators are bare. 

To write down explicitly the SC equation in the closed form, we must adopt the 
approximation or ansatz for the truncation of infinite hierarchy of SC equations in general 
field theories. Once such an ansatz (e.g., quenched ladder approximation in the SD 
equation) is adopted, we should not make further approximation to solve SC equations. 
However this forces us to solve the infinite set of integral equations distinguished by an 
integer n in ITF. In most studies based on the SC equation in IFT, therefore, infinite sum 
over discrete Matsubara frequencies was evaded by searching for frequency-independent 



approximate solution from the first |T9|, ^ |21|, |2^, To simplify further the 

equation, the constant (i.e., momentum-independent) solution has been investigated as 
in the NJL model in the leading 1/N approximation. By using such solutions, some 



qualitative features of FTD phase transition have been investigated in QED4 [|^, ^ 



QED3 [P, m, H and QCD4 H, 13 • However it is evident that in gauge theories the 



zero-temperature limit of the momentum-independent approximate solution obtained in 
such a way does not necessarily coincide with the zero-temperature solution [|l5l. To 



avoid this type of discrepancy, we must obtain the momentum-dependent solution at 
finite-temperature as tried in for QED3. 



On the other hand, Thermo Field Dynamics (TFD), a realization of Real-Time For- 



malism (RTF) of FTD field theory proposed by Takahashi and Umezawa ||28|, ^ 
needs no discrete sum and treats only continuous variables. Hence TFD is well suited to 
extend the usual field theoretical methods to the finite-temperature and is expected to be 
quite efficient to write down and solve the SC equation. Thanks to TFD, various methods 
developed so far in solving the integral equation can be extended to SC equations in the 
FTD case without much difficulty. Actually the infinite set of SC equations indexed by 
an integer in ITF is reduced to the corresponding single SC equation with an additional 
argument of continuous variable in TFD, as exemplified for the SD equation in this paper. 
Such a viewpoint has been overlooked to the best of our knowledge. 

This paper is the first of a series of papers which treat the transition of finite- 
temperature and density in field theories by use of self-consistent equations based on RTF, 
particularly on TFD. In this paper, restricting to the Nambu-Jona-Lasinio (NJL)[^-type 
four-fermion model and QED, we write down the SD equation for the fermion propagator 
based on TFD. Based on the solution of the SD equation for the fermion propagator, 
we show the existence of chiral-symmetry-restoring transition in QED4 by obtaining the 
critical line which separates the low-temperature phase where the chiral symmetry is spon- 
taneously broken from the high-temperature phase where the chiral symmetry restores. 
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This paper is organized as follows. In section 2, the formahsm of TFD is briefly 
reviewed to prepare the necessary materials and to fix the notation. In section 3, we treat 
the simplest case of NJL-type four-fermion model in the leading 1/N approximation, 
i.e., chain approximation. Within this approximation both formalism give the same SD 
equation for the fermion mass function. In section 4, we write down the SD equation in 
FTD QED based on TFD. In section 5, we solve the SD equation under the instantaneous- 
exchange approximation. In section 6, we consider another approximation and compare 
the result with the previous one. The final section is devoted to the conclusion and 
perspective. 
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2 Thermo-Field-Dynamics 

In this section we recall basic materials of TFD which are necessary to derive the SD 
equation. 



2.1 scalar 

Corresponding to the lagrangian of free scalar field (at zero temperature), 

>C = ^(5,0)^-^mg0^ (3) 
the free scalar propagator at zero-temperature reads 

In the TFD, the free scalar propagator at finite-temperature is given by the matrix form 



32, 33 



where 



withQ 



zD'^\k) = U{(3, k) (^1^^) J^^^ ) f/(/?, k), (3 ^ 1/T, (5) 

TT(a i\ f cosh 9k sinh6'fc\ 
^(^'^) = Uinh4 cosh^j' 

1 e-/3|feol/2 
cosh6'i. := , - — - , sinh^fc := - — - . (7) 

In the zero-temperature limit P oo, Uab{P, k) 6ab and hence the scalar particle A(k) 
and the thermal ghost A*{k) decouple. 

The explicit form of the propagator reads 

abf,, fcosh^ Ok A{k) +smh^ Ok A* {k) cosh^^ sinh4[A(A;) + A*(fc)] \ 
^ ' Uosh4sinh^fc[A(A;) + A*(A;)] sinh^ 9kA{k) + cosh^ 9kA*{k) J ' ^ ' 

Using the identities, cosh^ 'dk — sinh^ -(9^ = 1 and cosh^ -dk + sinh^ = 2 cosh^ — 1 
= coth^^^, we obtain 

iD^\k) = cosh^ OkAik) + sinh^ ekA*{k) 
= A{k) + sinh^ ek[A{k) + A*{k)] 
= A{k) + 2'K6{k^ - ml)NB{k) 

^ +2TT6{k'^ -ml)NB{k) 



k"^ — ttiq + ie 

+ 7r6{k'' -ml) coth^^, (9) 



-I- /I I t-"^ ^1 

k"^ — rn? 



"0 



-'^The notation A := B implies that A is defined by B. 
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where Nsik) is the Bose- Einstein distribution function: 



Nsik) := sinh^ 6, = / (10) 

exp{(3\ko\) - 1 



The full propagator will be assumed to be of the form 

zV^^k) = U{f3, k) (^(^) ^ ) [/(/?, kl (11) 

with T'(fc) being a complex function. 
As a consequence of the SD equation: 

tV''\k) = iD''\k) + iD^''{k){-iW^)iV'^\k), (12) 

the self-energy matrix is written as [| 

-^Ii'^\k) = u-\p,k)[^-'^^^^^ zn*V))^"'^^'^)' ^^2) 



where each component is related as follows |33 



ni2(fc) = n2i(A;) = -2icosh^feSinh^fcSn(fc) = -ztanh2^fcSn^^(A;), 

Ii^\k) = -Yi^^*{k). (14) 

This allows us to write the function T>{k) as 

Hence the SD equation for the scalar propagator reads 

iV-\k) = i/\-\k)-Yi{k), (16) 



where the real and imaginary part is related to 11^^ (p) as 

mi{k) = m^\k), 

<^U{k) = e{ko)tanh[^kopIl^\k). (17) 

2.2 fermion 

For the lagrangian of free fermion field, 

C = ^/j{t ^-mo)^, (18) 

^This is proved if there is the spectral representation which imphes the Kubo-Martin-Schwinger (KMS) 
equihbrium condition |p5[. 
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the free fermion propagator is obtained as 



S(P) = — ^ = r (19) 



Then the free fermion propagator at finite-temperature (and -density) is given by 



where S*{p) denotes 



= = -» r^V"^ . (21) 

p — niQ — le p^ — ttlq — le 



and 

cos -e(po)e~^^/^ sin yjp. 

, e (po ) e^''/^ sin v?p+^ cos 

with e(po) := 6'(po) - ^'(-po) and 



k(ap,/.) = ( ^"^^^"'^ ) , (22) 



_ g(po)e/^(p°+^)/4 + g(-po)e-/^(P°+^)/^ 
"^"^ "^"^^ ~ Ve/5(P"+A')/2 + e-/5(ft.+M)/2 ' ^^^^ 

_ 0(po)e"^(*'°+^)/^ + ^(-po)e^(P°+'')/^ 
^'"^ '^^+'^ ~ Ve/^(Po+M)/2 + e-/3(po+A.)/2 ■ ^^^^ 

In the zero-temperature hmit /? — > oo, V"afe(/?,p, jj) — > and hence the thermal part S*{j)) 
is separated. 

Exphcitly writing down the matrix element: 

iS''\p) 

cosVp,M'5'(p) - sinVp,M'S'*(p) -e(po)e"^''/^ cos V5p_^ sin v9p,^[S' + -S"*] (p).^^. 
e(pp)e/3/^/2 cos (/?p,^ sin ¥?p,^[S + S*]{p) - sin^ ¥5p,^^(p) + cos^ ¥?p,^5*(p) 7 

Note that, since cos^ ^p+ii + sin^ Vp+pL = 1, and cos^ fp+fi — sin^ V^p+^ = 2 cos^ V^p+/i — 1 = 
e(po) tanh ^^^^^y^, the 1-1 component of the fermion propagator is rewritten as 

iS^'^ {p) = cos^ (Pp+^lS (p) - sin^ '-Pp+^S* (p) 
= S{p)- sin^ ipp+f,[S{p) + S*{p)] 
= S{p) ~2ne{po)S{p^ -ml){fi + mo)NF{p) 

^ - ' ^ -'2 2 



2ne{po)6{p' - m^)(^ + mo)NFip) 



— rriQ + le 

+ 7re(po)5(p'-mg)(^+mo)tanh^^^^-t^, (26) 



^ — mo 

where Np{p) is the Fermi-Dirac distribution function: 

NM := sinVp+. = ^,(p„+t) ^ / (Po) + ^-,(p„^,) ^ ^ ^(-Po). (27) 



The full propagator will have the form 



with S{p) being a complex function. 

Similarly in the case of the scalar field, the SD equation: 

iS'^^p) = is^\p) + iS''^{p){-ij:"^)iS'^\p), 

is compatible with the self-energy matrix written as 
-zE^\p) = V-\P,p) ^-'fP^ 

where 

= -e-^^^J:'\p) = te{po)e-^^/He^n2^p+,QJ:'\p), 
Y?\p) = 

This allows us to write the function S{p) as 



S{p) 



IfiPf, - mo - S(p) + ie' 
Hence the SD equation for the fermion propagator takes the form: 

is-\p) = ts-\p)-j:{p), 

where 

m{p) = ^j:^\p), 

aS(p) = e{po)coth[^{po + ix)p^''{p). 

2.3 photon 

Given the lagrangian for the photon field 

with a being the gauge-fixing parameter, the photon propagator is given 

d 



-9tiu - (1 - 0!)kf,K 



where D {k)\jn=o is the massless scalar propagator. 



For the full photon propagator 'D^t(^) obeying the SD equation 
the vacuum polarization function is introduced 



(37) 



(38) 



where 



nil(A;) = Uf,{k) = -2i cosh 0k sinh 0k^U^.{k) = -i tanh 2^fe$5ni^,(A;), 



nS(A;) = -nii.*(A;). 
It is expected to satisfy the relation 

iV;;l(k)^iD-Xk)-u^,{k), 



(39) 



(40) 



where 



m^,{k) = mf,{k), 

^n^.{k) = e(A;o)tanh[^A;o]S>n;i(fc). 



(41) 



2.4 remark on the full propagator 



In this paper we discuss only the real part of the self-energy function as explained in 
section 4. 

In the general case of Q'S(p) 7^ 0, the full boson propagator should take the form 



F - n(A;) + ie 
and the full fermion propagator 



k'^-U{k)+ie k^-U*{k)- 



le 



Nsik), (42) 



^- E(p) + ie 



^— S(p)-|-ie ^— E*(p)— ie 



Nf{p). (43) 
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3 SD equation for NJL model 

We consider the model with four-fermion interaction of Nambu-Jona-Lasinio (NJL) type 
whose lagrangian is given by 

C = ^-^ + (a=l,...,iV) (44) 



which is equivalent to 



(45) 



where we have introduced the scalar and the pseudoscalar auxiliary fields a and vr respec- 
tively. 

Due to the Yukawa interaction of the fermion with the auxiliary scalar field, two 
tadpole diagrams contribute to the fermion self-energy in the NJL model in the leading 
order of 1/N expansion (Fig.l), which is the chain approximation: 



iii 



(2^ 



{D"(0)tr[z5^i(p)] + Di'(0)tr[icS22(p)]} 



(46) 



where D°^(}z) is the auxiliary field propagator. Since Z)^^(O) = G/N , D^^(O) = 0, it turns 
out that the second diagram has vanishing contribution. 
When mo = 0, we notice that 



ii[iS^\p)\ =itr(l) 
Hence we obtain 

where 



^ + z6(po)vr5(p2 - M2)M tanh ^i^i±i^ 



p2 _ ^2 



(47) 



-vr 



GtT{l)e{po)S{pl - El)M tanh 



P{Po + yu) 



-vrGtr(l) - : + M tanh + , (4J 



2Ef 



Ep := VWTW, P = \P\. 
Then the integration with respect to the time- component is straightforward 

dpo „trfl) M 



f ^D'\0)mr[^S'\p)] = -G- 



4 Ep 

Identifying the fermion mass M with the real part of S 0: 

M = tr(3?S)/tr(l) = tr(3?Sii)/tr(l) 



tanh ^^^^ + tanh ^i^^ 



(49) 



(50) 



(51) 



•^See the next section for the reason. 
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therefore, the gap equation in the NJL^ {D > 2) model at FTD is obtained 

^D-lp ^ 



(27r)^-i Ep 



2 2 



(52) 



The solution obtained from this gap equation is momentum-independent. 

This gap equation indeed coincides with that obtained by performing the sum over 
Matsubara frequencies in ITF. 



n=+oo „ iD—ip 
n=—oo V " / 



(53) 



where 



S{P) 

m 



1 





^ - E ~ p2 _ 5]2 

G 
N' 



, po = (2n + l)i7rr + /i, 



(54) 



The self-consistent solution E is p-independent and hence we can perform the discrete 
sum: 

oo ^D-lp 5^ 



E = tr(l)GT Y: J 



n=—oo 



„ jU-lp y 



tanh^i^^ + tanh^i^^ 



(55) 



with Ep := VP^ + E^. 
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4 QED at finite-temperature and -density 



It should be mentioned on treatment of imaginary part of self-energy function in this 
paper. In TFD, two types of field indexed by 1 and 2 appear in the theory where the 
type-1 field is the usual field and the type-2 newly introduced tilde field corresponding to 
the ghost field in the heat bath [0^ . 

The self-energy function differs from S^^(p), although the real parts of and 
S^^(p) coincide p3| : 



KS(p) = 3?S"(p), = e{po) coth (56) 

In this paper we consider only the real part of the self-energy 3fJS(p). In other words, we 
search for a self-consistent solution satisfying Q'E(p) = 0: 

^E{p) = 3?S^^(p), = 0. (57) 

This is one possible solution in the scheme of SC equations. This standpoint is different 
from the finite-temperature perturbation theory []33| , |3^ . An imaginary part of the self- 
energy at finite temperature describes the approach to equilibrium and can be related 
to the dissipative transport coefficients of viscosity and heat conductivity In this 



paper we consider only the equilibrium case and neglect the imaginary part which will be 
discussed in the subsequent paper. R 



4.1 derivation of the SD equation from TFD 

We consider the fermion self-energy function S^^(p) in QED in D-dimensions (QED/j). 
In the bare vertex approximation, the first non-trivial contribution to S^^(p) is written 
by using only S^^ and Vj^l as diagrammatically shown in Fig. 2: 



(2vr)^5(p - q - kh,tS''{qh,tVllik) 



{2n)D J {2tt) 



(5^ 



since the 1-1 component S^^ of the full fermion propagator corresponds to two physical 
type-1 external legs. 

In this paper, moreover, we search for the self-consistent (real) solution of the full 
fermion propagator in the following form: 



' - M{p) 



+ 7re{po)6{p^ - M\p)) tanh ^i^2±i^(^+ M{p)) 



(59) 



*The limit 9S(p) should be taken from the negative side SS(p) < 0. 
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and we choose the photon propagator in the Feynman gauge at finite-temperature: 



+ 7i6{k^ -U{k)) coth 



-g^,iV^\k), (60) 



where n(A;) is the vacuum polarization function of the photon. Then the real part of 
T}^{p) is given by 



{27t)d J (27r) 



(27r)"5(p - q - k)Z{q)^,{4+ M{q)h^ 



re{qo)6{q^ - M\q)) tanh 



5(A;2 -n(A;))coth^ 



F-n(fc) ' q^-M^{q) J" ^^^^ 

Here we define the wave function renormalization function Z{p) and the fermion mass 
function M{p) through 



Hence we obtain 



M{p) 



i iZ{p) 
~ ^- ~ ^ - M{p) ' 

tr[g?S(p)] _ tr[a^S"(p)] 
tr(l) ~ te(T) ' 
^ tr[^3?S(p)] 



(62) 



tr[^3fJS"(p)] 



p2tr(l) p2tr(l) 
Thus we can write down the coupled SD equation in FTD QED^ as 



(63) 



Z-\p) = l~{2~D)elix 



d'^qp-q^,, r e(go)5(g' - M2(g)) tanh[f (go + /i)] 



(27r)^ p2 



+ 



{p - qY - Yi{p - q) 
6{{p - qf - n(p - q)) coth[f Ipo - q 



q^-M'^{q) 



M{p) 



-Deln 



d^'q 
(27r)« 



Z{q)M{q) 



+ 



6(go)^(g2-M2(g))tanh[f(go + /i)] 
{p - qf - n(p - q) 

5((p-g)2-n(p-g))coth[f|po-go|] 



It is easy to see that if M{p) is a solution of the SD equation, then —M{p) is also a 
solution. 

In ITF we must solve infinite number of integral equations (indexed by an integer) 
for the function M„(P) with one argument, see Section 5.1. In the TFD, on the other 
hand, we have only to solve the single integral equation for the function M(po, P) with 
two arguments. 

^ At zero temperature, the bare vertex approximation and the free photon propagator in the Landau 
gauge leads to no wave function renormahzation of the fermion propagator, Z{p) = 1, jist which is 
consistent with the bare vertex approximation in Hght of the Ward identity, Z\ = Zi- This property is 
not preserved at finite temperature in the exact sense. However the deviation of Z(j)) from 1 is at most 
logarithmic in the momentum and does not substantially change the result in this paper. This will be 
discussed in detail in the subsequent paper. 
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4.2 finite-temperature and -density QED4 

In what follows, we put Z{p) = 1 as explained in the footnote. For convenience, we 
decompose the self-energy function I]{pq, P;T, fi), right- hand- side of the above equation, 
into two parts, J2f{pQ, P;T, fi) and T,ph{pQ, P;T). Thus the SD equation for the fermion 
mass function M{pq, P) is written as 

M(po, P) = S/(po, P; T, /x) + Sp,(po, P; T). (66) 

In what follows, we replace for simplicity the vacuum polarization function of the photon 
with its infrared value, i.e., mass of the photon generated by finite-temperature and - 
density effects p6[| . 

n(fc) =^ n(0) = m^ = nv'iT, ij,). (67) 

Improvement of this approximation will be tackled in the subsequent paper. 
The first self-energy part reads 

f d^'q , ,e(go)5(g'-M2(g))tanh^to) 
Ef(po,P;T,fi) = -eluD / —^M(qo,Q) ^ — . 68 

Decomposing the integration measure into the angular part and the radial one, we obtain 
for Q D > 3 

J^fipo, P; T,fi) = - cI-kDCd ^ r Q^'-^dQ d^ sin^-=^ ^ 

j~oo Ztc Jo Jo 

[qo, L^) _ ^^^2 _ (p2 + g2 + ^2) + 2PQ cost?' ^ ' 

where P = \P\, Q = \Q\, cosd = P-Q/PQ and 

2 7l~ 1 

= (2vr)^-ir(^) = 2^-2vr^/2r(^)' ^ 

Here we have introduced the ultraviolet (UV) momentum cutoff A. 
For D = 4, after performing the angular integration, we obtain 



nj r+00 pA 

S/(po,P;T,/i)= -/ dqo dQt{q,)6{q^ - M\q)) 

TT J-00 JO 



M(go, Q) tanh[-(go + /.)]- In _ _ _ _ , (71) 



where we have introduced a new coupling constant (C4 = l/(47r^)) 



= 3 case (QED3) will be discussed in a separate paper. 
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Next we carry out the integration with respect to go in the formal way. Assuming that 
the equation 

pl-P'-M'ipo,P)=0 (73) 
has two solutions: [] Po = i-E'p {E$ > 0), 

S{p'-M\po,P)) 

= 5{pI-P^-M\po,P)) 

5{po-Et>) , 5{po + Ep) 



Using this formula, the first self-energy part ^/(po, P] T, /i) reads 
^fiPo,P;T,fi) = - [ dQ^ 

IT Jo P 

M{E^,Q) t^nh^^^ ^Jpo-EQ^-iE^J' 



- WoM^i^o, Q)\,o=+E-^ iPo - EqY - (E^y 
M(-EQ,g) tanh^i^ ^Jpo + E^r - {E:^f 



2Eq - 4^M2(go, Q)\,,=.E- (Po + EqY - {E-f 



(75) 



where 



E^:=^\P±Q\^ + m^. (76) 
On the other hand, the second self-energy part Sp/i(po, P] T) reads {D > 3) 

2 , d^g (5((p-g)2-m2)coth^^^ 



Sp,(po,P;T) = -nDelJ—^M{qo,Q) 



{2n)D v-'-^ q2_Q2_M^{q,,Q) 

r,r 2 r^dq, f^^o-2,n ■ , M(go,Q)coth^i^ 

J-oo 27r Jo Jo % ~ Q ^ (O'o, Q) 

5{iPo - qof - {P^ + Q^ + m") + 2PQ cos^). (77) 

Note that T,ph{po, P',T) has no /i-dependence. 

For D = 4, after performing the angular integration, we obtain 

Eph{po,P;T) = / dQ dqo — — — ^ . 78 

TT Jo Jv P qQ - - M2(go, Q) 



Note that the range of integration — 1 < cos < 1 is equivalent toP = {J\P — Q]"^ + m? 



< \pq — go|< \ {P + QY + "^^} from the fact that the delta function has the support at 



^We assume that such a pair of solutions exists. Note that invariance of the mass function M{pQ, P) = 
Ml—po, P) under the reflection pa — po docs not hold in the presence of the chemical potential fJ. ^ 0. 



If the solution has no po-dependence, M {po, P) = M (P), then po = ±Ep, Ep J P^ + RP{P) 
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a point COST? = P^-^Q-'+^^HPo-q,? ^ r^^^^ obtain 



n Jo r 



-I 



dqoM{qo, Q)^ — -f- — 

dqoM{qo, Q)^ — -f- — 



(79) 



Integration with respect to qo can not be performed for 'Ephipo, P; T), unless we know the 
dependence of the solution M{qo,Q) on go- 
Thus the SD equation of FTD QED4 for the fermion mass function is written as 



TT Jo 



M{El Q) tanh (p, _ EqY - 



qo=+E+ 



M(-i?Q,g) tanh ^^(to) {po + EQf-{Et,f 

{Po + Egy - {E-f ^ 



coth 



+-/ dQ^ / dqoMiqo,Q)-, 
n Jo P [Jpo-E- %-Q - M^{qo, Q) 

-f 

Jpo 



rp^+E+ ^ ^ cothta^ 

dqoMiqo, Q)^ — -j- — 



(80) 



4.3 limiting cases of QED^ 



The SD equation in the various limiting cases is obtained in the following replacements. 
In the zero-temperature limit /5 — > 00, 

tanh ^t^^Zif) ^ 0^Eq - /x) - ^(/x - Eq), 
tanh ^t^^±i^ ^ 1 ^ e(EQ - /.) + ei^i - Eq), 



coth 



1. 



(81) 



In the zero-temperature /3 — > 00 and the zero-density limit — > 0, 

1, 



tanh^i^^,tanh^i^^ 



coth^^^-1. 



(82) 
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5 Instantaneous-exchange approximation 

5.1 instantaneous- exchange approximation in ITF 

In the imaginary-time formalism, the full fermion propagator as the solution of the SD 
equation has the form (when the wave function renormalization can be neglected): 



S~\p) =^ + S„(P), po = (2m + l)nT - tfi, 



(83) 



which becomes frequency-dependent. Then we can not perform the discrete sum in ad- 
vance, before we find the solution of the self-consistent SD equation: 



E^P) = elT J ^^-^D,,{ko =p,-qo,K = P-Q) 



S„(Q) 



^4) 



with go = (2^ + l)vrT — ifj,. This fact renders the analytical treatment as well as the 
numerical one extremely cumbersome. 

To avoid this problem, the instantaneous-exchange (IE) approximation p2[] is adopted: 

D^^{ko,K)-D^^,{ko = 0,K). (85) 

In this approximation the solution Sm(P) becomes frequency-independent S(P) and the 
summation over discrete frequencies can be performed explicitly as 



2(P) = el 



Q 



where 



^6) 



Eq = JQ' + I:'{Q). 



(87) 



For our choice 



the angular integration is performed for > 3 as 



S(P) = DCoel r g^-'dQ 

JO 



s(g) 



4E, 



Q 



X / M- 



sin^-3 d 



p2 + g2 + ^2_2pg cos^9' 
Then the SD equation for QED4 in IE approximation is obtained 



(88) 



a 



S(P) = - dQ 



TC Jo 



, gs(g) 

PEo 



tanh^^^ + tanh^i^ 



In 



{P + Qf + m^ 
{P-Qf + m?' 



^9) 



(90) 
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5.2 instantaneous- exchange approximation in TFD 

We consider the similar approximation of the SD equation in TFD, which we call also 
the IE approximation. If we put /cq = in the photon propagator, iVj^Kko, K), the 

temperature-dependence piece vanishes from the nature of the delta- function 5(— il'^ — 
m^) = as long as > 0, i.e.. 



tVlliko = 0,K) = -9,u^. (91) 

Therefore, in our IE approximation 

Ep;,(po,P;T) = 0. (92) 

Furthermore, Sj(po, P-, T,fi) becomes po-independent, so that the fermion mass function 
M{pq,P) obtained in the self-consistent way should be po-independent, M{P). Thus the 
SD equation in the IE approximation reads 

M P = — / dQ—- — tanh + tanh In— — ^ . 93 

In the IE approximation our choice of the photon propagator gives the same SD equation 
in the two formalism; ITF and RTF (or TFD). 

In the zero-temperature limit (at finite-density /z 7^ 0), the SD equation reads 

^'(^) - ^ £ <'«7i^M(Q)«(^« - i^!g:':r . (94) 

5.3 numerical results {fi = 0) 

We have solved the integral equation: 

^<^> = ^ r ^Q/I^^^W) ^^-''if (P ! g)^ + 3 ■ ^« = V<3^ + AW. (95) 

where the photon mass is borrowed from the 1-loop calculation ||36| : 



ml = le^T^ = ^aT\ (96) 

In the numerical calculation the dimensionful quantities e.g., T, P,Q, M{P), Eq, ■ ■ ■ 
is normalized by the cutoff A, and hence the SD equation can be rewritten in terms of 
dimensionless variables where the range of integration is the finite interval [0, 1]. 

First of all, we study the zero-temperature limit T = 0: 

^^^^ = ^ lo 1^ [p-Qr ^ \/^'+^(Q)'- (97) 
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The fermion mass is dynamically generated due to the nontrivial solution M{p) 7^ and 
the chiral symmetry is spontaneously broken in the strong coupling region a > Oc with 
the critical coupling ac = 0.39. 

Fig. 3 shows the coupling- const ant a dependence of the fermion mass M(0) which 
is identified with a pole position of the fermion propagator: p"^ = M'^{p'^) ~ M^(0). 
Fig. 4 exhibits the fermion mass function M{P) which is monotonically decreasing in P. 
These features are qualitatively the same as the well-known results obtained from the 
zero-temperature SD equation in QED [1^. 

At finite temperature T > 0, Fig.5 corresponds to Fig. 3 at zero-temperature. The 
critical coupling ac depends on the temperature and adT) is greater than ac at zero- 
temperature. Temperature-dependence of M(0) for a fixed a{> ac) is exhibited in 
Fig. 6 (a). This implies that there exists a critical temperature above which the chiral 
symmetry restores in the sense that the dynamical fermion mass vanishes: M(0) = for 
T>Tc. 

In Fig. 7 we compare the fermion mass functions M{P) for various temperature T at a 
fixed a. For relatively low-temperature, there exists a peak for M{P). However the peak 
disappears as the temperature increases and the critical temperature is approached. 

The phase diagram for two-coupling space (a,T) is depicted in Fig. 8 where the low- 
temperature strong coupling phase where the chiral symmetry is spontaneously broken is 
separated by the critical line from the high-temperature chiral-symmetry-restoring phase. 
This shows that no matter how large the coupling constant a is taken there exists a finite 
critical temperature < 00, namely, the chiral symmetry always restores at a finite 
temperature T^. 

To see the scaling of the dynamical fermion mass, the dimensionless quantity M(0)/Tc 
in the region T/Tc ~ 1 is plotted in Fig.9 for various values of 9 := {a — a^/ac, a = 
a(T ^ 0) = 0.36. Fig.9 suggests the existence of the scaling function, i.e., M{0)/Tc 
is written as a single function of T/Tc near the critical temperature Tc. Our data are 
still insufficient to specify the critical exponent associated with the finite-temperature 
transition. More detailed data will be given in the subsequent paper. 

Finally some remarks on the numerical calculation are in order. In our numerical 
calculation, we adopt the Double Exponential (DE) formula to choose sample points. DE 
formula is efficient when there is a weak singularity at the end of the interval of integration. 
Actually the integral kernel of the SD equation at T = has a singular point at P = Q. 
At finite-temperature T 7^ 0, on the other hand, the singularity at P = Q disappear due 
to the existence of the photon mass rriT- Therefore we can calculate the integrand also at 
P = Q without difficulty for high-temperature. In our numerical calculations we need at 
least 200 sample points according to DE formula to guarantees the precision of 10~^ for 
the solution M{P). However we need more sample points to obtain the stable result in 
the relatively low-temperature region. As the number of sample points becomes larger, 
the low-temperature result becomes more stable, i.e., the rapid change of M(0) at T ~ 
disappear even in the low-temperature region (T/A < 0.01) as shown in Fig.6(b). In order 
to check the influence of the singularity to a minumum, we divided the interval [0,1] into 
subintervals so that each of the subinterval contains 50 sample points which are chosen 
according to DE formula. For more details, see |3^ 
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6 po-independent approximation 



In this section we consider another approximation to the SD equation. In the IE approx- 
imation the second self-energy part Tiph including the factor coth[|go] has automatically 
dropped from the SD equation and the mass function M(po, P) gets po-independent. To 
go beyond the IE approximation, we must include the effect of the second self-energy part 
Epfe. However the integral equation with two arguments: 

M(po, P) = S/(po, P; T, 11) + S,,(po, P; T), (98) 

is rather difficult to solve even in the numerical way. Therefore we consider po-independent 
solution M(P) of the SD equation by requiring 

S^M)^,. (99) 

We try to determine the po-independent solution M{P) self-consistently so as to satisfy 
this condition. The previous IE approximation is nothing but a sufficient condition for 
the solution to be po-independent. 



6.1 finite-temperature case (fi = 0) 

Prom the g'o-independence of the solution, the first self-energy part reads 



Ef{po,P;T,0) = — / dQ-^M{Q) tanh 

ZTT Jo i^^Q 



X 



X 



In - - (g); + In + - 



with 



(Po - EqY - {E-y 



Eq - ^g2 + M(g)2. 



(po + EqY - {E-y 



(100) 



(101) 



Here the value of po is specified below. Note that the right-hand-side of the above equation 
can be rewritten in the form which is invariant under the refiection: po — > —po'- 



E;(po,P;r,0) 



a 



27r 
X In 



.A 

/ dQ 

Jo 



Q 



PEc 



M{Q) tanh '■ — - x 

\pl-{E^ + EQfU~{E^-EQf] 

\pl-{E- + EQy\\pl-{E--EQy] 

On the other hand, the second self-energy part reads 



(102) 



^ph{po,P;T) 



a 



27r Jo 



''dQ-f^M{Q)If,{P,Q,M], 



PE, 



Q 



(103) 
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where 



Ip{P,Q,M] 



Po-E„ 
PO-E+ 



+ 



Po+E„ 



dqo ^ 

dqo 



Qo + Eq go — Eq 
( 1 1 



coth 



qQ + Eq Qq - Eq^ 



coth 



(104) 



which is rewritten as 
UP,Q,M] -- 



dt 



+ 
2 



dt 
dqo 



t-po- Eq t-po + Ec 
1 1 



coth — 
2 



t + po + Eq t + po- Eq^ 

Qo + Eq go - Eq 



coth 



/3t 



{qo + EQY-pl (qo-EQ) 



Po 



coth 



(105) 



Hence T,ph{pQ, P;T) is also invariant under the replacement: Pq —>■ —po- In the zero- 
temperature limit (T = 0), 



Spft(po, P] 0) 



a 

Ytc Jo 



\q-R-M{Q) 



In 



[pl-iE^ + EQnpl-iE;n-EQf 

H-iE;;. + EQnpl-iE:^-EQy 



(106) 



In the zero-temperature and zero-density limit {T = = n), therefore, the SD equation 
M(P) = S^(po, P; 0, 0) + Sp;,(po, P; 0) reads 



M(P) 



TT Jo 



Q M{Q) ^^pI-{P + Q + E^ 



Q) 



P Eq pI-(\P-Q\+Eq) 



(107) 



Now we discuss how to choose po- To determine po, we require positivity of the 
nontrivial solution M{P) obtained self-consistently, at least in the neighborhood P 0, 
which is satisfied if the integrand is positive for P ^ 0. In the zero-temperature limit 
T = 0, this follows from 



Po<Eq + Q. 
In the finite-temperature case T > 0, we notice 



(108) 



Ef{po,P;T,0) 



a 

¥tt Jo 
X In 



^ dQ-R-M{Q) tanh ^^"^ 



1 + 



+ C)(P' 



X 



1 



PEq 2 

4{Eq + y/g^ + m')QP 
[{Eq + v/Q' + m^)2 - piyQ'2 + m' 
AjEQ - VQ^T^)QP 

[{Eq - VQ^T^y - pI]VWT^ 



(109) 



®For the po-independent M{P), the integration with respect to go can be done in principle, but the 
result can not be expressed by the elementary function. 



22 



where we have used 



= ^Q' + m^± , + 0{P% (110) 

In the total (self-consistent) mass function M{P) = + Sp/^, we assume the first self- 
energy part T,f{pQ,P]T,0) is the dominant part. Hence T,f{po,P;T,0) has the same 
signature as M(P), since it is not difficult to show po cannot be chosen such that both 
E/(po, P — 0;T, 0) and T,ph{po, P ~ 0; T) are simultaneously positive (or negative). Thus 
positivity requirement of both the total mass function M{P) and the first self-energy part 
T,f{pQ,P ~ 0; T, 0) restrict the range of po as 



0<EQ~JQ^ + m^<Po<EQ + JQ^ + m^. (Ill) 



Here Eq > \^Q'^ + m?^ i.e., M{Q) > m = ^-^aNfT is satisfied for sufficiently strong 
coupling a > ftc ^ 1 and sufficiently low-temperature T/A <C 1 in the case of finite fermion 
fiavor Nf, while it is always satisfied in the quenched limit Nf 0. This condition does 
not guarantee positivity of the second part Sp/i(po, P] T) even at T = 0. The apparently 
simplest choice po = is excluded by this condition. 

6.2 numerical results = 0) 

In the numerical calculation we take the simplest choice: 

Po = Eq, (112) 

and study the effect of the second part of the self-energy contribution. For this choice, the 
SD equation for the mass function reads M(P) = Sj(i?Q, P; T, 0) + Tiph^Eq, P; T) which 
can be rewritten in the form: 

M(P\ 2a Q r^t. [gg - 2^^] tanh[f Eq] - EqQq coth[f gp] 

M{P) = — / dQ——M[Q) / rfgo W^r-^^ ' ^3 

vr Jo PEq Je- go Wo - (2^q)^J 



where 



47r 



E^ := ^|P ± g|2 + ml, ml = —aT\ (114) 



At first, we study the zero-temperature limit. 



MlP\ " f^n <2 ,,„nu {P + Q)(P + Q + 2Eq) 

M(P) = dQ—M(Q) In (115) 

The chiral symmetry is spontaneously broken in the strong coupling region a > ac = 0.62 
due to dynamical generation of the fermion mass M(0), as shown in Fig. 10. The dynamical 

^The fact that S/(po: -P; T) and Tiph{po, P] T) has an opposite sign is consistent with the result of the 
perturbation theory at finite temperature p3|. 
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fermion mass function M{P) is shown in Fig. 11. These results are quahtatively in good 
agreement with those in the IE approximation. 

Next, we study the finite-temperature case. Fig. 12 shows the temperature-dependence 
of the fermion mass M(0). The critical temperature Tc exists at Tc ~ 0.0651, above which 
the chiral symmetry restores. 

Temperature-dependence of the fermion mass function is shown in Fig. 13. This shows 
the same tendency as the case under IE approximation. 

6.3 comparison with results in two approximations 

Finally we compare two results obtained in IE approximation and in po-independent 
approximation. 

Fig. 14 is the plot of the dynamical mass M(0) obtained at T=0 under two approxima- 
tions where the horizontal axis is the coupling constant normalized respectively by each 
critical coupling. This shows that M(0) under po-indep. approximation is slightly smaller 
than IE approximation for all over the range of the coupling constant. In other words, 
additional effects from Ep/j work to decrease M{P) at T=0. 

Fig. 15 shows the plot of the M(0)/Tc versus T/Tc for the same ratio a^T = 0)/a = 
0.618. Hence M (0)/Tc in IE approximation exhibits larger value for T/Tc ~ and smaller 
value for T /Tc — 1 than po-indep. approximation. 

Concerning the transition temperature, IE approximation gives slightly smaller critical 
temperature {Tc/A = 0.063) than po-indep. approximation {Tc/A = 0.065). This implies 
that effect of the retarded propagation of the photon has a tendency of raising the critical 
temperature. 
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7 Conclusion and Discussion 



As the formalism which can treat the field theory at finite-temperature and finite-density, 
we know two formalism: Imaginary-Time Formalism (ITF) and Real-Time Formalism 
(RTF). In this paper we have derived the Schwinger- Dyson (SD) equation at finite- 
temperature and non-zero density as a tool to analyze the non-perturbative effect in 
the Nambu-Jona-Lasinio (NJL) model of four-fermion interaction and QED. In section 3, 
we have shown that two formalism give the same SD equation for the NJL model in the 
leading order of 1/N expansion, i.e., chain approximation. 

In D-dimensional QED (QED^)) at finite-temperature and nonzero- density we have 
derived the SD equation for the fermion propagator under the bare vertex ansatz (ap- 
proximation). This is a coupled integral equation for the fermion mass function M(po, P) 
and the wave function renormalization function Z{pq, P) with two- variables po, P = \P\. 

In order to obtain the explicit solution, we have put Z{p) = 1 from the outset and 
searched for the real solution of the SD equation for the fermion mass function M{p) 
in QED4 at finite-temperature and nonzero-density. Even in this stage, solving the two- 
variable integral equation is rather difficult even in the numerical calculation. Therefore 
we consider the approximation so that the SD equation reduces to the one-variable inte- 
gral equation whose solution is denoted by M(P). Actually we have adopted two types 
of approximations: instantaneous-exchange (IE) approximation and the po-independent 
approximation. The IE approximation in TFD adapted in this paper leads to the same 
SD equation as that derived under the corresponding IE approximation in ITF In IE 
approximation only one temperature-dependent factor tanh[|i?p] appears in the SD gap 
equation and hence all the effects of temperature comes from this factor. Moreover IE 
approximation greatly simplifies the gap equation, since IE approximation automatically 
render the solution M{p) po-independent: M{P). 

The SD equation derived in this paper based on RTF immediately allows us to do an- 
other approximation, i.e., po-indep. approximation which goes beyond IE approximation. 
This approximation requires, from the first, that the solution M{p) is po-independent. 
However this approximation is able to incorporate a new temperature-effect coming from 
the factor coth[|go], which is disregarded in IE approximation. 

Numerical calculations under the two approximations have shown qualitatively the 
same result. For any coupling constant a{> ttc) belonging to the strong coupling region 
where the chiral symmetry is spontaneously broken in the zero-temperature, there exists 
a finite critical temperature Tc{a) < 00 above which the chiral symmetry restores. No 
matter how the coupling a may be strong, the critical temperature Tc is kept finite. 

In this paper we have not analyzed the finite-density transition. In the presence of 
chemical potential 7^ 0, the symmetry pq —pQ will be lost in the SD equation 
and we are forced to tackle with the original two-variable integral equation. Therefore 
approximations adopted in this paper must be improved to treat such a case. This will 
be a subject in the forthcoming paper. 
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9 Figure Captions 

Fig.l: Fermion self-energy contribution in the NJL model. 

Fig. 2: Fermion self-energy contribution in QED. 

Fig.3: M(0) vs 1/a at T = in IE approximation, = 0.39). 

Fig.4: Fermion mass function at T = in IE approximation, [a = 0.40). 

Fig.5: M(0)/A vs 1/a at T/A = 0.0794 in IE approximation, {a^ = 0.70). 

Fig. 6: (a) Temperature-dependence of the fermion mass in IE approximation (a — 0.63, 
Tc/A — 0.0628), (b) Dependence of M(0) on the number of sample points N in the 
low-temperature region T/A < 0.01. 

Fig. 7: Temperature-dependence of the fermion mass function in IE approximation for 
T/A = 10-^ 10-^ 10-^ 10-^ 0.04 {a = 0.70). 

Fig.8: Phase diagram of QED4 at finite-temperature in IE approximation. 

Fig.9: M{0)/Tc vs T/Tc in IE approximation where 6 = 0.08, 0.75, 0.94, 1.7. 

Fig. 10: M(0)/A vs 1/a at T = in po-indep. approximation, (ac — 0.618). 

Fig. 11: Fermion mass function at T = in pQ-indep. approximation, (a = 1.0). 

Fig. 12: Temperature-dependence of the fermion mass in po-indep. approximation {a — 1.0, 
Tc/A = 0.0651). 

Fig. 13: Temperature-dependence of the fermion mass function in po-indep. approximation 
for T/A = 10-^ 10-^ 0.024 {a = 1.0). 

Fig. 14: Comparison of IE and po-indep. approximations: M(0)/A vs ac/a. 

Fig. 15: Comparison of IE and po-indep. approximations: M(0)/Tc vs T/T^ for ac/a = 
0.618. 
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